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We report the first observation in a high energy neutrino telescope of cascades induced by at- 
mospheric electron neutrinos and by neutral current interactions of atmospheric neutrinos of all 
flavors. Using data recorded during the first year of operation of IceCube's DeepCore low energy 
extension, a sample of 1029 events is observed in 281 days of data. The number of observed cascades 
is A'cascadc ~ 496 ± 66(stat.) ± 88(syst.) and the rest of the sample consists of residual backgrounds 
due to atmospheric muons and charged current interactions of atmospheric muon neutrinos. The 
flux of the atmospheric electron neutrinos is determined in the energy range between approximately 
80 GeV and 6 TeV and is consistent with models of atmospheric neutrinos. 

PACS numbers: 95.55.Vj, 14.60. Lm, 29.40.Ka, 95.85.Ry, 25.30.Pt 



The spectrum of atmospheric Vfj, (we do not differenti- 
ate between 17 and ly here) has been measured at energies 
up to 400 TeV Much less is known about atmo- 
spheric Ve. These electron neutrinos come mostly from 
the decays of kaons and muons produced in cosmic-ray 
air showers. Underground water Cherenkov telescopes 
like IMB-3 and Super-Kamiokande as well as calorimet- 
ric detectors such as Frejus, NUSEX, and Soudan-2 have 
studied atmospheric with energies up to a few tens 
of GeV but no measurement has been made at 

higher energies. So far, searches for at higher ener- 
gies have yielded upper limits Theoretical flux 
calculations for atmospheric Ve are poorly constrained at 
energies above 100 GeV due to the uncertainties in kaon 



production ll|, [12 1 . 

Atmospheric neutrinos with energies below 100 GeV 
are of particular interest for studies of neutrino oscilla- 
tions. The first oscillation maximum for i/^ — i>r, and 
corresponding minimum in the v^, survival probability, 
occurs at 24 GeV for vertically upward-going neutrinos 
13[. More energetic atmospheric are an important 



no-oscillation baseline for comparison with lower ener- 
gies, and are also an irreducible background for charm- 
induced and astrophysical neutrino searches. 

In this Letter, we report on a measurement of atmo- 
spheric neutrino-induced cascades using the DeepCore 
infill array in IceCube. Here, the "cascades" refer to 
i>e charged current (CC) interactions and neutral cur- 
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FIG. 1. Schematic side view of IceCube. The DeepCore infiU 
array is shown in the shaded region. 

rent (NC) interactions of neutrinos of all flavors. From 
a selected sample of cascade candidate events, an atmo- 
spheric Ve flux is obtained in the energy range between 
about 80 GeV and 6 TeV. 

IceCube is a high-energy neutrino detector buried in 
the Antarctic ice. It observes Cherenkov light from neu- 
trino interactions. The main detection signatures for 
neutrinos are long, straight tracks and approximately 
spherical cascades. The former are created by neutrino- 
induced muons while the latter are produced by neutrino- 
induced electromagnetic and/or hadronic showers. The 
DeepCore inflll array [14|] to IceCube reduces the energy 
threshold of IceCube to energies as low as 10 GeV. Deep- 
Core's denser optical module spacing, higher quantum ef- 
ficiency photomultiplier tubes, and lower trigger thresh- 
old, along with its deployment in the ice with the best 
optical properties in the lower half of the IceCube detec- 
tor, all enhance low-energy neutrino detection. 

This analysis used the first IceCube data run with 
DeepCore, from June 2010 to May 2011 when a total of 
79 IceCube strings, including six specialized DeepCore 
strings, were operational. Each string consists of 60 dig- 
ital optical modules (DOMs), equipped with a photo- 
multipler tube [l^ and data acquisition electronics [l6| . 
The DeepCore fiducial volume included these six strings 
and the seven adjacent standard IceCube strings. Since 
then, the remaining seven IceCube strings have been de- 
ployed, including two additional DeepCore strings. The 
schematic view of IceCube is shown in Fig. [T] 

The DeepCore fiducial volume contains the 454 DOMs 
deployed at depths greater than 2100 m on the 13 strings 
of DeepCore. A dedicated DeepCore trigger was run on 
these DOMs. It read out the full IceCube detector if pho- 
tons ("hits") were observed in local coincidence (LC) on 
at least three neighboring DOMs within 2.5 /Lts. Details 
of the LC circuitry are available in Ref. [l^ . The average 



trigger rate was 185 Hz and was a factor 13 smaller than 
the total IceCube trigger rate. 

To avoid observational bias, 10% of the data, dis- 
tributed evenly through the year, were used to develop 
the analysis and verify the detector simulation. The re- 
sults presented here are based on the remaining 90% of 
the data set. After application of general data quality 
criteria, the 281 live-days of data that remain are used 
in this analysis. 

To observe cascades, one must reject two types of back- 
grounds which are more numerous than the desired sig- 
nal. The first class of background is muons produced in 
cosmic-ray air showers in the Earth's atmosphere ("at- 
mospheric muons"), which penetrate the 2100 m of ice 
above DeepCore and produce hits in the fiducial vol- 
ume. The second class, much less common than the first, 
consists of atmospheric also produced in air show- 
ers, which undergo CC interactions in the ice but which 
produce relatively low energy muons. Simulated data are 
obtained from Monte Carlo (MC) programs modeling the 
detector response to both of these types of background, 
as well as neutrino-induced cascades. An extensive air 
shower simulation 1 1711 i s used for atmospheric muons and 
a separate program [ISj is for neutrinos weighted with the 
Honda [n| and the Bartol [3| atmospheric flux predic- 
tions. Selection criteria are applied to the data to reduce 
these backgrounds sufficiently to observe cascades. Fig. [2] 
compares the performance of these criteria with the pre- 
dicted performance. 

Following the trigger, an online veto algorithm identi- 
fies events where at least one LC hit is observed in the 
IceCube veto region (i.e., all IceCube sensors outside the 
fiducial region). If the times of the veto hits relative 
to the hits in the fiducial volume are consistent with a 
relativistic particle traveling thr( 3ug h the veto region to 
DeepCore, the event is rejected [l^. This filter reduces 
the rate of atmospheric muons in DeepCore by a factor of 
11 with respect to the DeepCore trigger, with negligible 
losses of neutrino interactions within the fiducial volume. 
An additional trigger criterion is applied in software, con- 
sidering not only the LC hits but also non-LC hits which 
are not too isolated in time or space from other hits in 
the event. The trigger demands at least eight hits in 
the event, of which at least four must be in the fiducial 
volume. 

Next, five variables quantifying the event topology are 
calculated. These are the depth of the first hit in the 
event, the sphericity of the event (the ratio of the small- 
est eigenvalue to the sum of all eigenvalues, the analogue 
of the tensor of inertia obtained by treating each hit as 
a point mass), the fraction of the total number of pho- 
toelectrons which are recorded within the first 600 ns 
(i.e., a measure of how fast the event develops), a sim- 
ilar fraction calculation excluding the two earliest hits 
assuming they may be due to noise, and the number of 
hits occurring in the veto region regardless of their time 
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or location. These variables are used to train a boosted 
decision tree (BDT) [g^. A cut on this 5-variable BDT 
reduces the number of data events to 2.7x10^, a factor 
of 1660 reduction from the DeepCore trigger. The num- 
ber of predicted atmospheric cascades is 2.3x10'* at this 
stage. 

More computationally intensive event reconstruction 
algorithms are run over the surviving events, succes- 
sively reconstructing the events under the hypotheses 
that they are produced by atmospheric muons or by 
neutrino-induced cascades. These likelihood- based re- 
constructions take into account the details of Cherenkov 
light propagation in the ice [2lj. A second BDT is then 
trained to discriminate between atmospheric muons and 
cascades, using seven variables: the ratio of the best- 
fit likelihoods obtained from these reconstructions, the 
depth of the first hit, the horizontal distance of the first 
hit from the center of DeepCore, the fraction of the to- 
tal number of photoelectrons detected within the first 
300 ns, and variables measuring whether the hit pattern 
tends to drift across the detector during the event, indica- 
tive of a muon-like track through the detector rather than 
a spherically expanding light pattern which is a typical 
signature of a cascade. The numbers of events passing 
a cut on the second BDT are 3.8x10'' (data), 2.1 xlO'' 
(atm. fi), 1.7x10^ (atni. CC), and 9.6x 10'^(atm. cas- 
cades). 

The data surviving this second BDT cut primarily con- 
sist of neutrino events, but the majority are produced by 
CC interactions rather than cascades produced by Ve 



or i>n NC interactions. 



Most of these CC events come 



from neutrinos interacting within the DeepCore volume. 
The hadronic cascade produced at the neutrino interac- 
tion vertex may obscure the muon track emanating from 
the vertex unless the muon has sufficient energy to travel 
a considerable distance and is oriented such that it does 
not escape the detector volume undetected. 

To reduce this CC background, we apply several 
additional criteria. We require that the reconstructed 
neutrino interaction vertex of the event not be near the 
top or bottom of the fiducial volume, where a muon track 
could be missed if it escapes into the uninstrumented 
ice below the detector or the relatively dusty ice above 
DeepCore, and that the first hit in the event also fall in 
this volume. We further require that the best-fit cascade 
likelihood be both relatively good and also better than 
that obtained from the track fit, that enough DOMs were 
hit in the event that the comparison between these fits 
should have discrimination power, and finally we require 
a high value for a high signal efficiency from the second, 
7-variable BDT 
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A total of 1029 events pass all of these criteria. We esti- 
mate from simulation that about half of these events are 
residual backgrounds due to atmospheric muons (14%) 
and CC events (36%). The remainder are Ve and 
Vfj. NC events (inseparable from i/e)- The MC predicts 
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FIG. 2. (Color online) The number of data events (black 
filled squares) passing each selection criteria is shown with 
the MC predictions: atmospheric muon (magenta open cir- 
cles), atmospheric Vfj, CC (red open squares), and cascade 
signal (blue filled circles). From left to right on the horizontal 
axis: all triggered events, DeepCore triggered events, Deep- 
Core filtered events, 5-variable BDT cut (BDT5), 7-variable 
BDT cut(BDT7), and the final selection. Lines are to guide 
the eye. 



TABLE I. The number of events in 281 days are shown after 
application of all selection criteria. N°'"^ denotes the number 
of observed data events. An average of the Bartol and Honda 
event rates is also shown. The neutrino simulations have sta- 
tistical uncertainties of less than 2% while the atmospheric 
muon MC has a statistical uncertainty of 45%. 
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that half of the cascades are produced by (primarily 
CC interactions) and half are NC events. The final 
event rates are shown in Table HI For atmospheric muons, 
the overall rejection is 4x10* with respect to the total 
IceCube trigger rate. However, the limited size of the 
simulated atmospheric muon sample which could be pro- 
duced with the available computational resources results 
in a large statistical uncertainty in this background. 

The largest systematic uncertainty in the signal predic- 
tion comes from the light detection efficiency in a DOM 
in situ. Varying the efficiency by 10% in the simula- 
tions, the predicted atmospheric neutrino rate changes by 
11% for i^f^ and 10% for Vg- The same procedure is per- 
formed for the atmospheric muons, except at an earlier 
stage (BDT7) of the analysis to ensure adequate statistics 
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and gives 30% uncertainty. The systematic uncertainties 
due to the optical properties of ice are estimated as 8% 
for atmospheric muons and 6% (2%) for atmospheric i/^ 
(Uf.) by comparing final level rates from simulations with 
two different ice models. With in situ calibration light 
source data at all depths of the detector, the baseline ice 
model is derived from a global fit to the data obtaining 
a collection of parameters for scattering and absorption 
properties [131 , while the alternative uses fits to describe 
the time and amplitude distributions for emitter-receiver 
pairs |24i]. 

We conservatively estimate the uncertainty in the at- 
mospheric muon rate due to uncertainties in cosmic ray 
composition by comparing our baseline simulation, based 
on spectra of individual elements [i^ with a proton- 
only composition model. The comparison is made at the 
BDT7 stage to ensure sufficient statistics, and shows a 
rate variation of 25%. Additionally, a 20% uncertainty 
for the cosmic ray flux normalization and 6% for the sea- 
sonal rate variation are included. These are summed in 
quadrature and give a total of 33% cosmic-ray flux un- 
certainty. Though large, this cosmic-ray flux uncertainty 
is smaller than the statistical uncertainty in the atmo- 
spheric muon rate due to the limited MC sample, so we 
use this estimate as the systematic uncertainty. The sys- 
tematic uncertainty for neutrino-nucleon cross sections 
is estimated to be 6%. The atmospheric flux uncer- 
tainties of 9% are obtained by comparing the final event 
rates with the Honda and Bartol flux predictions. Neu- 
trino oscillations have a very small effect in this sample 
(1.8% for i/fj, and 0.1% for i^e) due to the relatively high 
energy, {WJ"") - 232 GeV (see Fig.^. The Vr contribu- 
tion is estimated to be less than 1% of the data sample 
assuming standard oscillation parameters [26j . The total 
systematic uncertainties are 14% (11%) for atmospheric 
i/^ (ve) and 45% for atmospheric muons, as shown in Ta- 
ble IE 

The atmospheric muon background and the atmo- 
spheric z/p CC background are subtracted from the 
data. The latter contribution is estimated by averaging 
the Bartol and Honda atmospheric neutrino predictions. 
Half the difference is included in the systematic uncer- 
tainties. We observe an excess of cascade events, 

iVcascadc = 496 ± 66(stat.) ± 88(syst.) , 

where the total statistical uncertainty includes statistical 
uncertainties of the two subtracted background compo- 
nents, and the total systematic uncertainty is a sum in 
quadrature of the CC systematic uncertainties and 
the atmospheric muon systematic uncertainties. Since 
part of the systematic uncertainties does not come from 
the flnal level comparison, we conservatively do not con- 
sider the correlations among the systematic uncertain- 
ties. The cascade signal has a signiflcance of 4.5 standard 
deviations. We estimate based on simulations that 240 
± 66(stat.) ± 109(syst.) of the cascades are produced 



TABLE II. Systematic uncertainties. 
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FIG. 3. (Color online) The event rate as a function of the re- 
constructed cascade energy. The sum of all MC expectations 
(green) is consistent with 281 days of data rate. The dotted 
lines show the Bartol prediction while the solid lines indicate 
Honda predictions for the atmospheric neutrinos. Data are 
shown with statistical errors. Systematic uncertainties (Ta- 
ble |II| are omitted for clarity. 

by i^e- The data appear consistent with Honda model, 
and slightly below Bartol model which predicts 127 more 
neutrino events in total. The lower rate prediction from 
the Honda model, especially in Ve, is due to the differ- 
ent treatment of kaon production in the atmosphere [l^] , 
and is shown in Table HI Both Honda and Bartol estimate 
roughly 15% uncertainties in the at mosp heric Ve flux at 
100 GeV rising to 25% at 1 TeV [III, H [ill ■ 

Likelihood reconstructions are performed on every 
event in the final sample, simultaneously fitting a cas- 
cade hypothesis for deposited energy and vertex position 
and time. A vertex resolution of 9 m and an energy 
resolution of 0.12 in logio(E/GeV) are obtained. The 
absolute energy scale uncertainty is found to be 0.1 in 
logio(E/GeV). Using the energy reconstruction in Fig. [3] 
(rebinned to get sufficient statistics and reasonable uncer- 
tainties in each bin), we subtract the atmospheric muon 
and the atmospheric z/^ CC and NC to estimate the 
i/g excess. The Ve excess is converted into fiux by normal- 
izing to the expected number of events from an average 
of the Bartol and Honda fiuxes. In each bin, the horizon- 
tal bar indicates the bin width. The marker is placed at 
a point representing the average reconstructed energy of 
the contributing events. The vertical error bars include 
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FIG. 4. (Color online) The electron neutrino spectrum (green 
open triangles). The conventional Ve (red line) and (blue 
line) from Honda, i/e (red dotted line) from Bartol and charm- 
induced neutrinos (magenta band) [281 ] are shown. Previous 
measurements from Super-K Frejus 0], AMANDA [13, 
[m and IceCube [J, [H are also shown. 

TABLE III. The E^$„ flux. is in GeV. 



logioE-'" 


- logioE-^^ 




Ej$,(GeV cm-^s-^sr-i) 


1.0 


- 2.2 


80 


(7.5 ±5.4) X 10"^ 


2.2 


- 2.6 


251 


(1.8 ± 1.4) X 10"^ 


2.6 


- 3.4 


865 


(4.1 ±3.1) X lO"*^ 


3.4 


- 4.6 


5753 





the statistical and systematic uncertainties (see Fig. [4| 
and Table [nj). 

In conclusion, we have observed atmospheric neutrino- 
induced cascades, produced by I'e CC interactions and 
NC interactions of all flavors in IceCube. The atmo- 
spheric Ve flux in the energy range between 80 GeV and 
6 TeV is consistent with current models of the atmo- 
spheric neutrino flux. Future work will reduce the sta- 
tistical and systematic uncertainties and allow the flux 
to be measured more precisely. The techniques used to 
identify neutrinos interacting within the DeepCore vol- 
ume and veto atmospheric muons originating outside the 
detector can be applied to studies of neutrino oscillations 
and searches for low energy neutrinos from astrophysical 
sources or from dark matter annihilation or decay. 
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